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Strong and electroweak corrections to the production of Higgs+2jets via weak 

interactions at the LHC 
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FIG. 1: Topologies for t-, u-, and s-channel contributions to 
qq ggH in LO, where q denotes any quark or antiquark and 
V stands for W and Z bosons. 
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Radiative corrections of strong and electroweak interactions are presented at next-to-leading order 
for the production of a Higgs boson plus two hard jets via weak interactions at the LHC. The 
calculation includes all weak-boson fusion and quark-antiquark annihilation diagrams as well as the 
corresponding interferences. The electroweak corrections, which are discussed here for the first time, 
reduce the cross sections by 5%, and thus are of the same order of magnitude as the QCD corrections. 
As argued in previous papers, where s-channel diagrams and interferences were neglected, the QCD 
corrections connected to interference effects are small. 
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INTRODUCTION 

The production of a standard Higgs boson in associ- 
ation with two hard jets is a cornerstone in the Higgs 
search both in the ATLAS [H and CMS [1] experiments 
at the LHC for the Higgs mass range between 100 and 
200 GeV, which is favoured by the global Standard Model 
fit to electroweak (EW) precision data. 

The production of Higgs-|-2jets receives two kinds of 
contributions at hadron colliders. The first type, where 
the Higgs boson couples to a weak boson that links two 
quark lines, is dominated by squared t- and u-channel- 
likc diagrams and known as the "vector-boson fusion" 
(VBF) channel. The hard jet pairs have a strong ten- 
dency to be forward-backward directed in contrast to 
other jet production mechanisms, offering a good back- 
ground suppression (transverse-momentum and rapidity 
cuts on jets, jet rapidity gap, central-jet veto, etc.). 
Upon applying ajDpropriate event selection criteria (see 
e.g. Refs. and Refs. 0,0 for more refer- 

ences) it is possible to sufficiently suppress background 
and to enhance the VBF channel over the second H-|-2jets 
mechanism that mainly proceeds via strong interactions. 
In this second channel the Higgs boson is radiated off 
a heavy-quark loop that couples to any parton of the 
incoming hadrons via gluons 3, 11 1 . According to a re- 
cent estimate [l^ hadronic production contributes about 
4—5% to the Higgs+2jcts events for a Higgs mass of 
120 GeV after applying VBF cuts. A next-to-leading or- 
der (NLO) analysis [ll| of this contribution shows that 
its residual scale dependence is still of the order of 35%. 

Higgs production in the VBF channel is a pure EW 
process in leading order (LO) involving only quark and 
antiquark parton distributions. Approximating the cross 
section by t- and w-channel diagrams only (without in- 
terference), because s-channel diagrams and interferences 
are rather suppressed, the corresponding NLO QCD cor- 
rections reduce to vertex corrections to the weak-boson- 
quark coupling. Explicit NLO QCD calculations in this 
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approximation [^, [l^ 14 , 1^ 1^ confirm the expectation 
that these QCD corrections are quite small, because they 
are shifted to the parton distribution functions (PDF) 
via QCD factorization to a large extent. The resulting 
QCD corrections are of the order of 5—10% and reduce 
the remaining factorization and renormalization scale de- 
pendence of the NLO cross section to a few per cent. 

In this paper we complete the previous NLO calcula- 
tions for the VBF channel in two respects. Firstly, we 
add the complete NLO EW corrections, and secondly we 
include all interferences in the QCD corrections. While 
all interferences are negligibly small, as expected, the EW 
corrections are of the same size as the QCD corrections 
and thus phcnomcnologically relevant. 



DETAILS OF THE NLO CALCULATION 

At LO, the production of Higgs-t-2jets via weak bosons 
receives contributions from the partonic processes qq 
Hgq, qq Hgg, qq Hgg. For each relevant configura- 
tion of external quark flavours one or two of the topolo- 
gies shown in Figure [T] contribute. All LO and one-loop 
NLO diagrams are related by crossing symmetry to the 
corresponding decay amplitude H qqqq- The QCD 
and EW NLO corrections to these decays were discussed 
in detail in Refs. ■ IT§|. in particular a representative 
set of Feynman diagrams can be found there. 

Evaluating 2^3 particle processes at the NLO level is 
non-trivial, both in the analytical and numerical parts of 



2 



the calculation. In order to ensure the correctness of our 
results we have evaluated each ingredient twice, result- 
ing in two completely independent computer codes yield- 
ing results in mutual agreement. The phase-space in- 
tegration is performed using multi-channel Monte Carlo 
techniques [l9| implemented in different ways in the two 
different generators. 



Virtual corrections 

The virtual corrections modify the partonic processes 
that are already present at LO; there are about 200 EW 
one-loop diagrams per tree diagram in each flavour chan- 
nel. At NLO these corrections are induced by self-energy, 
vertex, box (4-point), and pentagon (5-point) diagrams. 
The calculation of the one-loop diagrams has been per- 
formed in the conventional 't Hooft-Feynman gauge and 
in the background-field formalism using the conventions 
of Refs. ^] and [21j, respectively. The masses of the ex- 
ternal fermions have been neglected whenever possible, 
i.e. everywhere but in the mass-singular logarithms. 

In the s-channel diagrams intermediate W and 
Z bosons can become resonant, corresponding to WH/ZH 
production with subsequent gauge-boson decay. In order 
to consistently include these resonances, we implement 
the finite widths of the gauge bosons in the "complex- 
mass scheme", which was introduced in Ref. [22] for 
LO calculations and generalized to the one-loop level in 
Ref. [2^. In this approach the W- and Z-boson masses 
are consistently considered as complex quantities, defined 
as the locations of the propagator poles in the complex 
plane. The scheme fully respects all relations that follow 
from gauge invariance. 

The amplitudes have been generated with FeynArts, 
using the two independent versions 1 [lil and 3 [IBl. 
The algebraic evaluation has been performed in two com- 
pletely independent ways. One calculation is based on 
the in-house Mathematica program that was already used 
in the algebraic reduction of NLO corrections to the 
H 4 fermions decays 17, [3]. The other has been 
completed with the help of FormCalc [2^. 

The tensor integrals are evaluated as in the calcula- 
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tion of the corrections to e+e~ 4 fermions [23|, [27 1. 
They are recursively reduced to master integrals at the 
numerical level. The scalar master integrals are evalu- 
ated for cornplex masses using the methods and results 
of Refs. m, [2^, [s^l • Tensor and scalar 5-point functions 
are directly expressed in terms of 4-point integrals [31 1. 
Tensor 4-point and 3-point integrals are reduced to scalar 
integrals with the Passarino-Veltman algorithm [32] as 
long as no small Gram determinant appears in the reduc- 
tion. If small Gram determinants occur, the alternative 
schemes described in Ref. 33| are applied. 



Real corrections 

Real QCD corrections consist of gluon emission and 
processes with gq and gq initial states. Analogously real 
photonic corrections comprise photon bremsstrahlung 
and photon-induced processes with 7(7 and 75 initial 
states. The matrix elements for these corrections have 
been evaluated using the Weyl-van der Waerden spinor 
technique as formulated in Ref. [33 | and have been 
checked against results obtained with Madgraph [sst . 

All types of real corrections involve singularities from 
coUinear initial-state splittings which are regularized 
with small quark masses. The mass singularities are ab- 
sorbed via factorization by the usual PDF redefinition 
both for the QCD and photonic corrections (see, e.g., 
Ref. ^^). Technically, the soft and collinear singulari- 
ties for real gluon or photon emission are isolated both 
in the dipole subtraction method following Ref. 37| and 
in the phase-space slicing method. For gluons or photons 
in the initial state the subtraction and slicing variants de- 
scribed in Ref. {36l | are applied. The results presented in 
the following are obtained with the subtraction method, 
which numerically performs better. 



NUMERICAL RESULTS 



We use the input parameters as given in Ref. [17|. 
Since quark-mixing effects are suppressed, we set the 
CKM matrix to the unit matrix. The electromagnetic 
coupling is fixed in the scheme, i.e. it is set to 
ac^ = V^GfiMy^^sl^/n, because this accounts for electro- 
magnetic running effects and some universal corrections 
of the p parameter. 

We use the MRST2004QED PDF ^] which consis- 
tently include 0{a) QED corrections. These PDF in- 
clude a photon distribution function for the proton and 
thus allow to take into account photon-induced partonic 
processes. As explained in Ref. [s^, the consistent use 
of these PDF requires the MS factorization scheme for 
the QCD, but the DIS scheme for the QED corrections; 
the corresponding factorization scales are identified with 
the Higgs mass Mh if not stated otherwise. We only 
use four quark flavours for the initial partons, i.e. we do 
not take into account the contribution of bottom quarks, 
which is suppressed. Since no associated LO version of 
the MRST2004QED PDF exists, we use these PDF both 
for LO and NLO predictions. For the renormalization 
scale of the strong coupling constant by default we em- 
ploy Mh, include 5 flavours in the two- loop running, and 
fix as{Mz) = 0.1187. 

Apart from the total cross section without any phase- 
space cuts, we consider the integrated cross section de- 
fined after applying typical VBF cuts to the outgoing 
jets. In this case, jets are defined from partons using 
the fcx-algorithm [s^, 0, 41 1 as described in Ref. [i^ . 
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More precisely, jets result from partons of pseudorapid- 
ity \v\ < 5 using the jet resolution parameter D = 0.8. 
We also recombine real photons with partons or jets ac- 
cording to this algorithm. Thus, some of the photons end 
up in jets, others are left as identifiable photons. Follow- 
ing Ref. [l5|, we specifiy the VBF cuts as follows. We 
require at least two hard jets with 



PTj>20GeV, |yj|<4.5. 



(1) 



where pxj is the transverse momentum of the jet and y-j its 
rapidity. The tagging jets ji andj2 are then defined as the 
two jets passing the cuts ([Ij with highest pt and pxji > 
PTj2- Finally, we demand a large rapidity separation of 
the two tagging jets by 



(2) 



In Figure [5] we plot the total cross section with and 
without cuts as a function of Mh. In the upper panel we 
show the absolute predictions in LO and in NLO includ- 
ing QCD and EW corrections. The VBF cuts reduce the 
cross section by a factor 3-4. In the lower panel we show 
the relative corrections. Without cuts the QCD correc- 
tions are about -1-5% and the EW corrections about —5% 
both depending only weakly on Mh and cancelling each 
other substantially. With cuts the EW corrections are 
approximately —6%, while the QCD corrections vary be- 
tween -3% and +2%. In the EW corrections the WW 
and ZZ thresholds are clearly visible. It is interesting to 
note that the EW corrections to the full VBF channel 
are similar in size and sign to the EW corrections to the 
subreactions pp WH/ZH -I- X Compared to the 
related decays the H WW/ZZ ^ 4/ [B, [i3| the size 
is similar, but the sign is different. 

In Table [T] we present integrated cross sections for 
Mh = 120, 150, 170, and 200 GeV without any cuts and 
in Table Ull results for VBF cuts. We list the LO cross 
section ctlo, the cross section ctnlo including QCD-I-EW 
corrections, and the relative QCD and EW corrections, 
Sqcd and Sew, respectively. The complete EW correc- 
tions Sew also comprise the corrections from photon- 



induced processes iJ-y-induced) which turn out to be > 
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and to reduce the EW corrections. We note that the 
QCD corrections are dominated by the known (vertex- 
like) corrections to the squared t- and u-channel VBF 
diagrams, while corrections to interference terms are at 
the level of 0.1%. These interference corrections are not 
enhanced by contributions of two t- or u-channel vector 
bosons with small virtuality and therefore even further 
suppressed when applying VBF cuts. 

In Figure [3] we show the dependence of the total cross 
section on the factorization and renormalization scale for 
Mh = 120 GeV. We set the factorization scale fi = /^f, 
which applies to both QCD and QED contributions, 
equal to the renormalization scale /ur = fj, and vary it 
between Afn/S and SAfn- In this set-up, we show the 
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FIG. 2: Higgs mass dependence of LO and complete NLO 
cross section (upper) and relative EW and QCD corrections 
(lower) without and with VBF cuts. 



Mh [GcV] 


120 


150 


170 


200 


O-LO [fb] 


5936(1) 


4271(2) 


3536(1) 


2743(1) 


(TNLO [fb] 


5890(2) 


4219(2) 


3538(1) 


2775(1) 


<5qcd [%] 


4.04(3) 


3.47(2) 


3.72(2) 


4.48(2) 


few [%] 


-4.81(2) 


-4.70(2) 


-3.65(1) 


-3.33(1) 


^-y-induced [%] 


0.86(1) 


1.04(1) 


1.14(1) 


1.27(1) 



TABLE I: Cross section for pp ^ H -I- 2jets -f X in LO and 
NLO without cuts and relative QCD and EW corrections. 



The contribution S. 



'7-inducC' 

is part of Sew) is also given separately. 



(J from 7-induced processes (which 



Mh [GeV] 


120 


150 


170 


200 


O-LO [fb] 


1830.5(5) 


1524.2(4) 


1353.8(3) 


1139.1(3) 


0"NLO [fb] 


1678.7(9) 


1422.9(7) 


1293.4(6) 


1106.0(5) 


Sqcd [%] 


-1.97(4) 


-0.60(4) 


0.41(4) 


1.76(3) 


few [%] 


-6.32(2) 


-6.02(2) 


-4.87(1) 


-4.64(1) 


'^T-induccd [%] 


1.14(1) 


1.21(1) 


1.25(1) 


1.31(1) 



TABLE II: As in Table U but with VBF cuts appUed. 
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FIG. 3: Scale dependence of LO and NLO cross sections with 
QCD or QCD+EW corrections for Mh = 120 GeV without 
cuts (upper) and with VBF cuts (lower), /xr = fj,F = fi for 
LO, NLO and NLO QCD, but hr = M^/^j, for NLO QCD'. 



LO cross section, the QCD corrected NLO cross section 
and the complete NLO cross section involving both QCD 
and EW corrections. In addition we depict the QCD cor- 
rected NLO cross section for the setup where /Ur = / fi 
(NLO QCD'). Varying the scale up and down by a fac- 
tor 2 (8) changes the cross section by 11% (29%) in LO 
and 3% (18%) in NLO for the set-up with VBF cuts. 
Without cuts the scale uncertainty it is at the level of 
3% (11%) in NLO, while it is accidentally small in LO 
for this specific Higgs-boson mass. 



CONCLUSIONS 

Radiative corrections of strong and electroweak inter- 
actions have been discussed at next-to-leading order for 
Higgs production via vector-boson fusion at the LHC. 
The electroweak corrections, which have not been cal- 



culated before, reduce the cross section by 5%, and are 
thus as important as the QCD corrections in this channel. 
QCD corrections to interference contributions turn out to 
be negligible, confirming previous approximations. 
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